#/usr/bin/Rscript

on_100_10 <- read.table("on_100_10/profits.dat")
#op_100_10 <- read.table("op_100_10/profits.dat")
#op_100_10 <- read.table("op_partial_fac5_100_10/profits.dat")
op_100_10 <- read.table("op_partial_fac5_100_10/profits.dat")
ov_100_10 <- read.table("ov_100_10/profits.dat")
ut_100_10 <- read.table("ut_100_10/profits.dat")
print("100 10 on")
shapiro.test(on_100_10$V1)
print("100 op")
shapiro.test(op_100_10$V1)
print("100 ov")
shapiro.test(ov_100_10$V1)
print("100 ut")
shapiro.test(ut_100_10$V1)

on_100_5 <- read.table("on_100_5/profits.dat")
#op_100_5 <- read.table("op_100_5/profits.dat")
#op_100_5 <- read.table("op_partial_fac5_100_5/profits.dat")
op_100_5 <- read.table("op_partial_fac5_100_5/profits.dat")
ov_100_5 <- read.table("ov_100_5/profits.dat")
ut_100_5 <- read.table("ut_100_5/profits.dat")
print("100 5 on")
shapiro.test(on_100_5$V1)
print("100 op")
shapiro.test(op_100_5$V1)
print("100 ov")
shapiro.test(ov_100_5$V1)
print("100 ut")
shapiro.test(ut_100_5$V1)


on_100_15 <- read.table("on_100_15/profits.dat")
#op_100_15 <- read.table("op_100_15/profits.dat")
#op_100_15 <- read.table("op_partial_fac5_100_15/profits.dat")
op_100_15 <- read.table("op_partial_fac5_100_15/profits.dat")
ov_100_15 <- read.table("ov_100_15/profits.dat")
ut_100_15 <- read.table("ut_100_15/profits.dat")
print("100 15 on")
shapiro.test(on_100_15$V1)
print("100 op")
shapiro.test(op_100_15$V1)
print("100 ov")
shapiro.test(ov_100_15$V1)
print("100 ut")
shapiro.test(ut_100_15$V1)


on_10_10 <- read.table("on_10_10/profits.dat")
#op_10_10 <- read.table("op_10_10/profits.dat")
#op_10_10 <- read.table("op_partial_fac5_10_10/profits.dat")
op_10_10 <- read.table("op_partial_fac5_10_10/profits.dat")
ov_10_10 <- read.table("ov_10_10/profits.dat")
ut_10_10 <- read.table("ut_10_10/profits.dat")
print("10 10 on")
shapiro.test(on_10_10$V1)
print("10 op")
shapiro.test(op_10_10$V1)
print("10 ov")
shapiro.test(ov_10_10$V1)
print("10 ut")
shapiro.test(ut_10_10$V1)

on_10_5 <- read.table("on_10_5/profits.dat")
#op_10_5 <- read.table("op_10_5/profits.dat")
#op_10_5 <- read.table("op_partial_fac5_10_5/profits.dat")
op_10_5 <- read.table("op_partial_fac5_10_5/profits.dat")
ov_10_5 <- read.table("ov_10_5/profits.dat")
ut_10_5 <- read.table("ut_10_5/profits.dat")
print("10 5 on")
shapiro.test(on_10_5$V1)
print("10 op")
shapiro.test(op_10_5$V1)
print("10 ov")
shapiro.test(ov_10_5$V1)
print("10 ut")
shapiro.test(ut_10_5$V1)

on_10_15 <- read.table("on_10_15/profits.dat")
#op_10_15 <- read.table("op_10_15/profits.dat")
#op_10_15 <- read.table("op_partial_fac5_10_15/profits.dat")
op_10_15 <- read.table("op_partial_fac5_10_15/profits.dat")
ov_10_15 <- read.table("ov_10_15/profits.dat")
ut_10_15 <- read.table("ut_10_15/profits.dat")
print("10 15 on")
shapiro.test(on_10_15$V1)
print("10 op")
shapiro.test(op_10_15$V1)
print("10 ov")
shapiro.test(ov_10_15$V1)
print("10 ut")
shapiro.test(ut_10_15$V1)

on_50_10 <- read.table("on_50_10/profits.dat")
#op_50_10 <- read.table("op_50_10/profits.dat")
#op_50_10 <- read.table("op_partial_fac5_50_10/profits.dat")
op_50_10 <- read.table("op_partial_fac5_50_10/profits.dat")
ov_50_10 <- read.table("ov_50_10/profits.dat")
ut_50_10 <- read.table("ut_50_10/profits.dat")
print("50 10 on")
shapiro.test(on_50_10$V1)
print("10 op")
shapiro.test(op_50_10$V1)
print("10 ov")
shapiro.test(ov_50_10$V1)
print("10 ut")
shapiro.test(ut_50_10$V1)

on_50_5 <- read.table("on_50_5/profits.dat")
#op_50_5 <- read.table("op_50_5/profits.dat")
#op_50_5 <- read.table("op_partial_fac5_50_5/profits.dat")
op_50_5 <- read.table("op_partial_fac5_50_5/profits.dat")
ov_50_5 <- read.table("ov_50_5/profits.dat")
ut_50_5 <- read.table("ut_50_5/profits.dat")
print("50 5 on")
shapiro.test(on_50_5$V1)
print("50 op")
shapiro.test(op_50_5$V1)
print("50 ov")
shapiro.test(ov_50_5$V1)
print("50 ut")
shapiro.test(ut_50_5$V1)

on_50_15 <- read.table("on_50_15/profits.dat")
#op_50_15 <- read.table("op_50_15/profits.dat")
#op_50_15 <- read.table("op_partial_fac5_50_15/profits.dat")
op_50_15 <- read.table("op_partial_fac5_50_15/profits.dat")
ov_50_15 <- read.table("ov_50_15/profits.dat")
ut_50_15 <- read.table("ut_50_15/profits.dat")
print("50 15 on")
shapiro.test(on_50_15$V1)
print("50 op")
shapiro.test(op_50_15$V1)
print("50 ov")
shapiro.test(ov_50_15$V1)
print("50 ut")
shapiro.test(ut_50_15$V1)

#QQplots

jpeg("normalidade_100.jpg")
par(mfrow=c(3, 4))

qqnorm(on_100_10$V1, main="on 100 10")
qqline(on_100_10$V1)
qqnorm(op_100_10$V1, main="op 100 10")
qqline(op_100_10$V1)
qqnorm(ov_100_10$V1, main="ov 100 10")
qqline(ov_100_10$V1)
qqnorm(ut_100_10$V1, main="ut 100 10")
qqline(ut_100_10$V1)

qqnorm(on_100_15$V1, main="on 100 15")
qqline(on_100_15$V1)
qqnorm(op_100_15$V1, main="op 100 15")
qqline(op_100_15$V1)
qqnorm(ov_100_15$V1, main="ov 100 15")
qqline(ov_100_15$V1)
qqnorm(ut_100_15$V1, main="ut 100 15")
qqline(ut_100_15$V1)

qqnorm(on_100_5$V1, main="on 100 5")
qqline(on_100_5$V1)
qqnorm(op_100_5$V1, main="op 100 5")
qqline(op_100_5$V1)
qqnorm(ov_100_5$V1, main="ov 100 5")
qqline(ov_100_5$V1)
qqnorm(ut_100_5$V1, main="ut 100 5")
qqline(ut_100_5$V1)

jpeg("homocedasticidade_100.jpg")
par(mfrow=c(3, 4))

plot(on_100_10$V1, main="on 100 10")
plot(op_100_10$V1, main="op 100 10")
plot(ov_100_10$V1, main="ov 100 10")
plot(ut_100_10$V1, main="ut 100 10")

plot(on_100_15$V1, main="on 100 15")
plot(op_100_15$V1, main="op 100 15")
plot(ov_100_15$V1, main="ov 100 15")
plot(ut_100_15$V1, main="ut 100 15")

plot(on_100_5$V1, main="on 100 5")
plot(op_100_5$V1, main="op 100 5")
plot(ov_100_5$V1, main="ov 100 5")
plot(ut_100_5$V1, main="ut 100 5")

jpeg("normalidade_10.jpg")
par(mfrow=c(3, 4))

qqnorm(on_10_10$V1, main="on 10 10")
qqline(on_10_10$V1)
qqnorm(op_10_10$V1, main="op 10 10")
qqline(op_10_10$V1)
qqnorm(ov_10_10$V1, main="ov 10 10")
qqline(ov_10_10$V1)
qqnorm(ut_10_10$V1, main="ut 10 10")
qqline(ut_10_10$V1)

qqnorm(on_10_15$V1, main="on 10 15")
qqline(on_10_15$V1)
qqnorm(op_10_15$V1, main="op 10 15")
qqline(op_10_15$V1)
qqnorm(ov_10_15$V1, main="ov 10 15")
qqline(ov_10_15$V1)
qqnorm(ut_10_15$V1, main="ut 10 15")
qqline(ut_10_15$V1)

qqnorm(on_10_5$V1, main="on 10 5")
qqline(on_10_5$V1)
qqnorm(op_10_5$V1, main="op 10 5")
qqline(op_10_5$V1)
qqnorm(ov_10_5$V1, main="ov 10 5")
qqline(ov_10_5$V1)
qqnorm(ut_10_5$V1, main="ut 10 5")
qqline(ut_10_5$V1)

jpeg("homocedasticidade_10.jpg")
par(mfrow=c(3, 4))

plot(on_10_10$V1, main="on 10 10")
plot(op_10_10$V1, main="op 10 10")
plot(ov_10_10$V1, main="ov 10 10")
plot(ut_10_10$V1, main="ut 10 10")

plot(on_10_15$V1, main="on 10 15")
plot(op_10_15$V1, main="op 10 15")
plot(ov_10_15$V1, main="ov 10 15")
plot(ut_10_15$V1, main="ut 10 15")

plot(on_10_5$V1, main="on 10 5")
plot(op_10_5$V1, main="op 10 5")
plot(ov_10_5$V1, main="ov 10 5")
plot(ut_10_5$V1, main="ut 10 5")

jpeg("normalidade_50.jpg")
par(mfrow=c(3, 4))

qqnorm(on_50_10$V1, main="on 50 10")
qqline(on_50_10$V1)
qqnorm(op_50_10$V1, main="op 50 10")
qqline(op_50_10$V1)
qqnorm(ov_50_10$V1, main="ov 50 10")
qqline(ov_50_10$V1)
qqnorm(ut_50_10$V1, main="ut 50 10")
qqline(ut_50_10$V1)

qqnorm(on_50_15$V1, main="on 50 15")
qqline(on_50_15$V1)
qqnorm(op_50_15$V1, main="op 50 15")
qqline(op_50_15$V1)
qqnorm(ov_50_15$V1, main="ov 50 15")
qqline(ov_50_15$V1)
qqnorm(ut_50_15$V1, main="ut 50 15")
qqline(ut_50_15$V1)

qqnorm(on_50_5$V1, main="on 50 5")
qqline(on_50_5$V1)
qqnorm(op_50_5$V1, main="op 50 5")
qqline(op_50_5$V1)
qqnorm(ov_50_5$V1, main="ov 50 5")
qqline(ov_50_5$V1)
qqnorm(ut_50_5$V1, main="ut 50 5")
qqline(ut_50_5$V1)

jpeg("homocedasticidade_50.jpg")
par(mfrow=c(3, 4))

plot(on_50_10$V1, main="on 50 10")
plot(op_50_10$V1, main="op 50 10")
plot(ov_50_10$V1, main="ov 50 10")
plot(ut_50_10$V1, main="ut 50 10")

plot(on_50_15$V1, main="on 50 15")
plot(op_50_15$V1, main="op 50 15")
plot(ov_50_15$V1, main="ov 50 15")
plot(ut_50_15$V1, main="ut 50 15")

plot(on_50_5$V1, main="on 50 5")
plot(op_50_5$V1, main="op 50 5")
plot(ov_50_5$V1, main="ov 50 5")
plot(ut_50_5$V1, main="ut 50 5")

dev.off()




